The role of static disorder in negative thermal expansion in ReOs 
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Time-of-flight neutron powder diffraction and specific lieat measurements were used to study the 
nature of thermal expansion in rhenium trioxide, an electrically conducting oxide with cubic sym- 
metry. The temperature evolution of the lattice parameters show that ReOs can exhibit negative 
thermal expansion at low temperatures and that the transition from negative to positive ther- 
mal expansion depends on sample preparation; the single crystal sample demonstrated the highest 
transition temperature, 300 K, and largest negative value for the coefficient of thermal expansion, 
a = — 1T(1) X 10~® K ~^ . For the oxygen atoms, the atomic displacement parameters are strongly 
anisotropic even at 15 K, indicative of a large contribution of static disorder to the displacement 
parameters. Further inspection of the temperature evolution of the oxygen displacement parame- 
ters for different samples reveals that the static disorder contribution is greater for the samples with 
diminished NTE behavior. In addition, specific heat measurements show that ReOa lacks the low 
energy Einstein-type modes seen in other negative thermal expansion oxides such as ZrW2 08. 



I. INTRODUCTION 



Metals that contract upon heating are exceedingly 
rare, if they exist at all. This observation can be 
understood by considering the types of cohesive in- 
teractions bonding atoms into solids. As a mate- 
rial is heated, interatomic distances expand due to 
the asymmetric shape of the potential energy well 
and therefore to anharmonic effects. [TJ |5] A chem- 
ical bond with a deeper potential well will confine 
atomic vibrations to remain more harmonic and con- 
sequently, covalently bonded materials such as di- 
amond exhibit lower thermal expansion than met- 
als such as copper. In some covalent crystals with 
open framework structures, the low thermal expan- 
sion can even become negative. Among the nega- 
tive thermal expansion (NTE) materials most stud- 
ied are framework oxides such as ZrW208[3l IH |5] 
and more recently metal cyanide frameworks such as 
Zn(CN)2,[51 17] all of which are electrical insulators 
or semiconductors. In contrast, the oxide ReOs is 
a superb electrical conductor with a room tempera- 
ture resistivity close to that of copper. [51 [S] Recently, 
Chatterji et al. have demonstrated through neutron 
powder diffraction and first-principles calculations 
that ReOa also exhibits isotropic NTE,|10j making 
it the only known simple binary oxide to combine 
metallic conductivity with NTE. 
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A topological feature that can lead to NTE behav- 
ior in an open framework structure is a linear chain 
consisting of a ligand coordinated to two metals or 
vice versa. [21 |TT] For example, simple binary oxides 
such as CU2O, Ag20 (0-M-O chains) [T^ and certain 
phases of Si02 (Si-O-Si chains) 1131114] demonstrate 
NTE behavior. While this is not a requisite for 
NTE behavior (magnetostriction and phase transi- 
tions can also result in NTE) , it is a simple principle 
by which new NTE materials can be pursued. The 
structure of rhenium trioxide contains such linear 
chains in a three dimensional framework consisting 
of corner-sharing ReOe octahedra (See Fig.[T|). The 




FIG. 1: The crystal structure of the cubic perovskite 
ReOs. ReOe octahedra connect via corner-sharing in 
infinite chains in all three directions. The atomic dis- 
placement ellipsoids of the oxygen atoms (in red) are 
represented at 90 % probability; their fiat shape indi- 
cates transverse displacement normal to the Re-O-Re 
bonds and into the void in the middle of the cubic cell. 
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simple ReOs structure is a prototype for closely re- 
lated materials such as the perovskite ferroelectrics, 
ferromagnets, and high Tc superconductors and can 
also be thought of as the cubic perovskite ABO3 
with a void in the place of the A cation. 

Thermal expansion investigations of ReOa have 
presented conflicting results. A laser interferome- 
try study on a single crystal concluded that ReOa 
displayed NTE between 100 K and 340K,[Tn]. In 
contrast to the single crystal work, a catalogue of 
powder X-ray diffraction (XRD) studies reported 
the coefficient of thermal expansion of ReOa to be 
small but positive for all temperature ranges. P^B] 
The latter study has been more frequently cited in 
reviews of NTE in oxide frameworks, [2], [17] leav- 
ing ReOa largely unrecognized as a potential NTE 
material until the recent article by Chatterji et al, 
which showed through diffraction studies NTE up to 
200 K.[lO^ In the present study, we show that ReOs 
can exhibit isotropic NTE up to room temperature, 
that the positive to negative transition is dependent 
on sample quality, and that specific heat measure- 
ments show a lack of any low energy Einstein-type 
modes that lead to the NTE behavior in other ox- 
ides. Furthermore, we show through careful analy- 
sis of the atomic displacement parameters the role 
of static disorder of the oxygen atoms in impeding 
NTE behavior and therefore explain the discrepancy 
between past thermal expansion studies of ReOa. 



II. EXPERIMENTAL METHODS 

The thermal expansion properties of various ReOa 
powder samples were investigated by performing 
variable temperature neutron powder diffraction 
(NPD) with time-of-ffight neutron data. The first 
sample investigated, labeled ReOs-a throughout this 
paper, was a powder prepared by the decomposition 
of an Re207-1,4 dioxane adduct at 140 °C as first de- 
scribed by Nechamkin et aZ. [18] The second sample 
studied, ReOa-b, was a powder purchased from Alfa 
Aesar (99.99% chemical purity). Finally, the third 
sample, ReOs-c, was prepared by crushing single 
crystals of ReOa grown by chemical vapor transport 
(CVT). By using HgCl2 as the transport agent in 
an evacuated Pyrex glass ampoule with ReOa pow- 
der (ReOa-b as the starting powder), the crystals 
were grown across a temperature gradient of about 
100 °C as described by Feller et aZ. [T3]. In addition 
to growing single crystals, CVT also separates the 
lower oxide phase impurities from the freshly grown 
ReOa crystals. Hence, sample ReOs-c was expected 
to be of higher chemical and phase purity than the 
others. 

NPD patterns were obtained on the high intensity 



powder diffractometer (HIPD) at the Lujan Neutron 
Scattering Center at the Los Alamos Neutron Sci- 
ence Center (LANSCE). HIPD is suitable for study- 
ing the behavior of materials as a function of tem- 
perature and pressure due to its relatively high data 
acquisition rates. HIPD data were collected on all 
samples between 15 K and 300 K. The Lujan Cen- 
ter employs a pulsed spallation neutron source, so 
the NPD histograms were collected in time-of-ffight 
mode on fixed detector banks. The backscattering 
detector banks, located ±153° with respect to the 
incident beam, provide the highest resolution pat- 
tern and higher-index reflections while the normal 
(±90°) and forward scattering (±40°) banks have 
a lower resolution but probe higher d-spacings. In 
the study by Chatterji et al, the neutron diffraction 
measurements were performed on a constant wave- 
length source (A = 1.359 ±0.001 A),\W^ limiting the 
diffractogram up to 8.5 A^^ in momentum transfer 
Q = 47rsin6'/A. More accurate structural param- 
eters such as atomic displacement parameters can 
be accessed with time-of-ffight neutron data since 
higher Q ranges can be probed. In our measure- 
ments, a Qmax of 16 was used for all the struc- 
tural refinements. 

The specific heat data were collected upon warm- 
ing from 2 K to 20 K on a 28.5 mg single crystal using 
a quasiadiabatic method as implemented on a Quan- 
tum Design Physical Property Measurement System 
(PPMS). Like sample ReOa-c, the single crystal was 
grown by CVT. 



III. RESULTS AND DISCUSSION 

Least squares structural refinements by the Ri- 
etveld method were performed with all six his- 
tograms of the NPD data using the GSAS software 
package. [20] In space group Pm3m, only four vari- 
ables were refined for each temperature: the lat- 
tice parameter; the isotropic displacement param- 
eter Uiso of Re, and the anisotropic displacement 
parameters Un (= U22), and U33 of O. To determine 
the lattice parameters of the sample accurately, X- 
ray powder diffraction was performed at room tem- 
perature with Si Standard Reference Material 640c 
purchased from National Institute of Standards (See 
Fig.|2]i). The neutron powder profile for the sample 
prepared from single crystals, ReOa-c, is shown in 
Fig.|2jD. In this sample, the lattice parameter plotted 
as a function of temperature clearly demonstrates 
NTE behavior from 15 K up to 300 K (See Fig.|3];). 
NTE in the other powder samples was indeed ob- 
served but limited to a narrower temperature range 
(See Fig.[3}i,b). For sample ReOa-a, the NTE behav- 
ior was observed only up to llOK, and for ReOa-b, 
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only up to 220 K. Clearly sample quality has an ef- 
fect on the thermal expansion properties and this 
may, in part, explain the discrepancy between past 
findings on the thermal expansion behavior of ReOs 
and the reason why Chattcrji et al. found NTE only 
up to 200 K from their neutron measurements while 
their dynamical lattice calculations predicted NTE 
up to 350K410^ 

The linear coefficient of thermal expansion for a 
cubic system a is defined as a = (1/a) x da/dT 
where a is the cubic cell parameter and T the tem- 
perature. A least squares polynomial fit (second- 
order) to the lattice parameters obtained from the 
neutron data affords a(T), which can be differenti- 
ated with respect to T and divided by the exper- 
imental values of a to obtain a as a function of 
T. For sample ReOa-c the average a below 300 K 
was found to be —1.1(1) x 10~^K~^. For sample 
ReOa-a, a was found to be -3.5(6) x lO-'^K"!, and 
-6.4(5) X lO-'^K-i for sample ReOs-b. Thus, low- 
ering of the transition temperature Tj^ite^pte is 
correlated with a smaller, negative value for a. 

For sample ReOs-c, the values for a{T) are close 
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FIG. 2: (a) Observed profile of ReOa from X-rays shown 
in open circles and the calculated profile as a smooth 
red line. The X-ray powder patterns contain both sam- 
ple ReOs-c and a Si standard. Upper tickmarks corre- 
spond to Re03-c and lower to the Si standard, (b) Ob- 
served profile from neutrons is shown in open circles and 
the calculated as a smooth red line for sample ReOa-c 
only with corresponding Bragg reflection tickmarks be- 
low. The difference curve between the observed and cal- 
culated patterns are represented by the smooth purple 
line below the powder profiles. 



to those determined by Matsuno et al. in their inter- 
ferometry study on single crystals where the value of 
a is —2 X 10^^ K^^ at about 120 K and slowly rises 
to -1 X lO"'^ at about 230K.[T5; While we did 
not make measurements above room temperature, 
Matsuno et al. found a to be positive and smaller 
than 2 X 10-6 K-i between 340 K and 500 K. [IF In 
an XRD study carried out on single crystals as well 
as powder samples, Chang et al. found a value of 
1.1(1) X IQ-^ K"^ above room temperature. [21] The 
volumetric a calculated by Chatterji et al. is not lin- 
ear as in our simplified model from the fitted curves, 
but nevertheless shows a volumetric a (= 1/3 linear 
a) varying between —2 x 10^® to zero below 350 
K.[Iir 

The atomic displacement parameters obtained 
from the NPD data offer some clues as to the mech- 
anism responsible for NTE in ReOa. Due to the 
high range of momentum transfer Q available from 
time-of-flight neutrons as well as the lack of any de- 
cay with Q in the atomic form factors for neutron 
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FIG. 3: Temperature evolution of the lattice parameters 
obtained from neutron diffraction for samples ReOs-a, 
ReOa-b, and ReOs-c. The smooth black curves are the 
least-squares, polynomial fits to the lattice parameters. 
Error bars are obtained from standard uncertainties in 
the structural refinements. 



4 



diffraction, accurate atomic displacement parame- 
ters were obtained from the Rietveld analysis (See 
Fig.|4^). The Re atom is located on the special po- 
sition la (site symmetry 3to), which constrains the 
atomic displacement parameter to remain isotropic. 
The oxygen atoms are located on the position 3c 
(site symmetry A/mmm), constraining them to have 
only two values, Uu (= U22) and C/33. While the dis- 
placement of the O atoms along the Re-O-Re direc- 
tion remains fairly constant, the transverse motion 
increases upon heating, as also observed by Chat- 
terji et aZ.[T^ The most striking feature is the gap 
between Uu (= U22) and C/33 as represented by the 
flatness of the displacement eUipsoids of the O atoms 
(See Fig.[T]). The anisotropy only increases with tem- 
perature, suggesting that the transverse motion con- 
tracts the cubic lattice upon heating and is therefore 
the most likely mechanism for the observed NTE be- 
havior. 

Further inspection of the oyxgen displacement pa- 
rameters of all three samples reveals differences that 
may explain the varying thermal expansion behav- 
ior. Though sometimes loosely referred to as ther- 
mal parameters or Debye- Waller factors, the values 
for Uij have contributions from both dynamic and 
static displacements of the atoms. Extrapolation of 
the Uij parameters for Re03 down to OK reveals 
non-zero values, which indicates significant contri- 
bution from static displacement, especially in the 
case of Uu (See Fig. [4^). The large gap between 
Uu and J733 at low temperatures was observed for 
all three samples. A useful parameter to describe 
the extent of anisotropy in the atomic displacements 
is the ratio of the maximum value to the minimum 
value of Uij]^^ in our case, it is U11/U33. When 
plotted versus temperature (See Fig.|4|3,c), this mea- 
sure of anisotropy reveals that it is constant (within 
the error bars) and large for sample Re03-a, while 
it is increases upon heating for Re03-b and Re03- 
c. Most importantly, the higher U11/U33 ratios for 
Re03-a and Re03-b than for ReOs-c at low temper- 
atures suggest more static disorder in the first two, 
which should affect the lattice dynamical properties. 
The M3 mode calculated by Chatterji et al involves 
rotation of the ReOg octahedra in the a — b plane 
and arises from transverse motion of the O atoms in 
this plane. [TU] Since the M3 mode has the largest 
contribution to the negative value for a at low tem- 
peratures, static displacement of the O atoms, trans- 
versely from the Re-O-Re bond, would diminish the 
NTE behavior. 

Specific heat measurements on a single crystal of 
Re03 were carried out to understand the type of 
lattice vibrations contributing to the specific heat. 
The low temperature behavior in the heat capac- 
ity of Re03 is that of a straightforward Debye sohd, 
with a Debye temperature O of 344 K and a Sommer- 



feld coefficient 7 of 2.8mJmol~^ (See Fig.|5^). 
These values agree very well with the values of 327 K 
for Q and 2.85 mJ mol"^ K^^ for 7 found in the early 
studies by King et al. 9 Fig.jsJ) shows the specific 
heat over vs. InT, which is an approximate rep- 
resentation of the one-dimensional phonon density of 
states of a solid. [251 A signature of an Einstein 
mode in this representation is a Gaussian-like peak- 
clearly missing in the low temperature region of the 
specific heat. These measurements strengthen the 
argument by Chatterji et al. from their lattice dy- 
namical calculations that NTE in ReOs is caused by 
low energy acoustic modes with the M3 mode con- 
tributing the most to the NTE behavior. [TO] Also, 
the specific heat results rule out the suggestion by 
Matsuno et al. that a low energy, bending-like opti- 
cal mode is responsible for NTE in ReO-i.[r5] 
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FIG. 4: (a) Temperature evolution of the atomic dis- 
placement parameters for ReOa-c, where Uiso repre- 
sents the rhenium isotropic displacement parameters and 
Uii,U22,U33 the oxygen anisotropic displacement pa- 
rameters, (b) The ratio U11/U33 for sample ReOa-a and 
(c) the same ratio for samples ReOa-b and ReOs-c. Uu 
and U22 represent atomic displacement of the oxygen 
atoms normal to the Re-O-Re bond, while f/33 to dis- 
placement along the bond. Error bars, shown within the 
symbols in (a), are obtained from the standard uncer- 
tainties in the structural refinements. 
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FIG. 5: Specific heat data for a single crystal ReOs. (a) 
Data fit below 4 K to the function C/T = "i + fiT^ (where 
= 12/57r'*i?//3; R is the gas constant). The fitting 
allows the Sommerfeld coefficient 7 = 3.6mJmol~^ 
and a Debye temperature of = 344 K to be obtained, 
(b) Data plotted as C/T^ vs. InT in order to reveal 
any low energy phonon spectral weight. The absence 
of any peaks indicates a lack of contribution from any 
Einstein-type modes; the Debye contribution appears as 
a constant up to and the rise at low In T is due to the 
electronic contribution (yT) of the delocalized electrons. 



Although both ZrW208 and ReOa have cubic 
symmetry, it is interesting to compare how their 
structures lead do different NTE behavior. The 
primitive cubic structure of ReOa (given the ab- 
sence of the A cation in its structure and thus an 
extremal value of the perovskite tolerance factor) is 
of course the consequence of metallicity: the large 
volume is maintained by the Fermi pressure of the 
delocalized d-electrons in the tt* conduction band. 
Nevertheless, the proximity to a 'buckling transi- 
tion' is revealed by a pressure-induced collapse of the 
structure produced by cooperative counter-rotation 
of the octahedra,[25j |26] in this case very similar 
to ZrW208.E71lll] The linear a obtained for ReOa 
is smaller than that observed in ZrW208, where 
a = —8.7 X 10~^ K~^. Indeed, the smaller ab- 
solute value of the NTE and the smaller tempera- 
ture range at which it occurs in ReOa compared to 
ZrW208 is consistent with the much weaker precur- 
sor effects in the pressure-induced collapse in the for- 
mer than in the latter. [25] In turn, the weakness of 
the fluctuations in ReOa is possibly a consequence of 



it possessing fewer soft degrees of freedom than the 
more open framework structure of ZrW208. In ad- 
dition, the low energy Einstein-type modes observed 
in ZrW208 and implicated in that material's NTE 
behavior, [SJ H3] are clearly missing in ReOa . 

The relatively small NTE effect observed in ReOs 
remains, nevertheless, remarkable because it is 
isotropic and unprecedented for a material with 
the electronic properties of a simple metal. While 
classifying an oxide as a simple metal is atypi- 
cal, the free-electron model works well enough in 
ReOs that it qualifies under the scheme of elec- 
trical transport properties (other examples include 
the Na^jWOa bronzes and Ru02).|29] Most simple 
metals have a much larger a than that of ReOa. 
For example, for copper, a is 16.5 x 10~^K~^ at 
room temperature. [30] Another notable exception 
is the intermetallic YGaGe, which has such small 
a that its behavior has been termed zero thermal 
expansion. [31] Like YGaGe, ReOa is also remarkable 
in that its a above room temperature is as low as 
in a highly covalent material like diamond where a 
is 1 X 10-s K-\[2l While ReOs does exhibit highly 
metallic behavior, the strong covalency between the 
Re d-states and O p-states and its unique frame- 
work topology allow for negative thermal expansion, 
demonstrating the propensity for transition metal 
oxides to challenge our conventions on bonding and 
properties in solid materials. Furthermore, ReOs's 
relatively simple structure has provided a way to 
quantify the role of static disorder of the oxygen 
atoms on interesting behavior such as negative ther- 
mal expansion. 
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